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A recently completed next-to-leading-order program to calculate neutrino cross sections, includ- 
ing power-suppressed mass correction terms, has been applied to evaluate the Paschos-Wolfenstein 
relation, in order to quantitatively assess the validity and significance of the NuTeV anomaly. In 
particular, we study the shift of sin^ 9w obtained in calculations with a new generation of PDF 
sets that allow s{x) s{x), enabled by recent neutrino dimuon data from CCFR and NuTeV, as 
compared to the previous s = s parton distribution functions like CTEQ6M. The extracted value of 
sin^ 9w is closely correlated with the strangeness asymmetry momentum integral x[s{x) — s{x)]dx . 
We also consider isospin violating effects that have recently been explored by the MRST group. The 
results of our study suggest that the new dimuon data, the Weinberg angle measurement, and other 
data sets used in global QCD parton structure analysis can all be consistent within the Standard 
Model. 
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Introduction: An important open question in parti- 
cle physics in recent years has been the significance of 
the "NuTeV anomaly" — a 3 a deviation of the measure- 
ment of sin^ 9w (0.2277 ± 0.0013 ± 0.0009) reported in 
Ref. P], from the world average of other measurements 
(0.2227 ± 0.0004). Possible sources of the NuTeV 
anomaly, both within and beyond the standard model, 
have been examined in j^]. No consistent picture has 
yet emerged in spite of extensive literature Jl, 0, IE 111 
on this subject. The measurement in Ref. |l| was based 
on a correlated fit to the ratios of charged and neutral 
current (CC & NC) interactions in sign-selected neutrino 
and anti-neutrino scattering events on a (primarily) iron 
target at Fermilab. This procedure is closely related to 
measuring the Paschos-Wolfenstein ratio 9], which pro- 
vides the theoretical underpinning of the analysis. Specif- 
ically, the Paschos-Wolfenstein ratio R~ is related to the 
Weinberg angle 9w by 

"cc "^cc 
~ - - sin^ 0w + '5i?;4 + ^Rqcd^ ^^EW 

where the three correction terms are due to the non- 
isoscalarity of the target (Ji?^), next-to-leading-order 
(NLO) and nonperturbative QCD effects {SRqqi^), and 
higher-order electroweak effects ((5i?^^). Since R^ is a 



ratio of differences of cross sections, the correction terms 
are expected to be rather small. But at the accuracy 
required to test the consistency of the Standard Model, 
all the corrections need to be quantified as precisely as 
possible — similar to some previous combined perturba- 
tive/nonperturbative re- analyses [3 El of challenges in 
QCD. 

In this paper, we focus on QCD corrections, which are 
generally recognized 0, 0, IE S to be the least well 
known. Let us write 

SRqcd = + SRJ + SR],Lo (2) 

where the three terms on the right-hand side are due 
to possible strangeness asymmetry (s^ — s ~ s ^ 0) 
and isospin violation (up^„ ^ dn.p) effects in the parton 
structure of the nucleon, and NLO (0(as)) corrections, 
respectively [2^. The original NuTeV analysis was car- 
ried out at LO in QCD and assumed SR~ = = SR7. 
Our analysis is based the recent NLO calculation of [12| , 
together with new parton analyses that explicitly allow 
strangeness asymmetry {6Rj 7^ 0) fl3| and isospin viola- 
tion {SRJ 0) jl^] . The actual calculation is carried out 
at the cross section level, i.e., using the first line of Eq. ^ 
rather than using the schematic linearized form given in 
the second line of Eq. and Eq. |5J|. Our results provide 
more realistic estimates of the sizes and uncertainties of 
the QCD corrections, and a new look at the significance 
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of the "anomaly." (Cf. also a recent re-evaluation of the 
electroweak correction to the calculation of R~ [l5|.) 

NLO Calculation: At sufficiently high neutrino en- 
ergy, the total neutrino cross section 
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can be calculated in QCD perturbation theory — in con- 
trast to charged lepton scattering, where the massless 
photon propagator leads to dominance of nonperturba- 
tive photoproduction events over deep inelastic scatter- 
ing. The differential cross section in Eq. © factorizes 
into a sum of convolutions of parton distribution func- 
tions and partonic cross sections 
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This calculation has been performed at NLO accuracy 
in Ref . ^3 • The analysis included target and charm mass 
effects. These corrections are needed to obtain reliable re- 
sults because there are non-negligible contributions from 
low Q values to the integral in Eq. Q — e.g., about 5% 
from < 1 GeV^ for crgc and from < 2 GeV^ for 
ctqq. (For NC neutrino events, it is not possible to ex- 
clude the low-Q region by experimental kinematic cuts.) 
Other corrections included are the non-isoscalarity of the 
target material (iron), i.e., SR'^ in energy averaging 
over the neutrino and anti-neutrino flux spectra; and cuts 
in hadronic energy (20 GeV < yE^, < 180 GeV for lepton 
inelasticity y) as used in the experimental analysis [l|. 

Ref. used previously available parton distributions 
[id. Il7j , all of which assume isospin symmetry and s — s 
symmetry within the nuclcon. The study confirmed the 
smallness|26j of the higher order corrections to i?^ in 
general. (The same conclusion is reached by the NLO and 
NNLO moment analyses of |EIEE|-) It was also shown 
that the non-monochromatic neutrino and anti-neutrino 
beams, with different profiles, and typical cuts in the 
hadronic event energy do not alter SRJ^j^q substantially. 
In the next two sections, we will examine shifts of the 
NLO calculation due to recent advances in global QCD 
analysis of parton distributions that allow strangeness 
asymmetry and isospin violation. 

In principle, the parton distribution functions in 
Eq. Q should be those of nuclear targets. Our calcu- 
lation is done as an incoherent sum of contributions from 
parton densities of unbound nucleons. This approxima- 
tion is reasonable in that we only calculate relative shifts 
between [S*"] = and [S'] ^ PDFs (where [S'] is de- 
fined in Eq. Q; similarly for isospin. In fact, experimen- 
tal information on nuclear PDFs is relatively scarce, and 
nuclear PDFs only account for leading twist 2 (r = 2) ef- 
fects. Higher twists, whether they relate to nuclear mod- 
ifications or not, are generally difficult to handle consis- 
tently. By limiting ourselves to r = 2, our error estimates 
may be underestimates. 



Strangeness Asymmetry: Because the strange quark 
mass rus is comparable to Aqcd , the strange quark PDF 
is a nonperturbative component of the nucleon bound 
state. Except for the strangeness number sum rule. 



[s{x) ~ s{x)] dx — , 



(5) 



there is no fundamental or approximate symmetry that 
relates the strange quark PDF s{x) to the antiquark PDF 
s{x). Limits on s~ = s{x) — s{x) can, therefore, only be 
derived from data (or perhaps eventually from a lattice 
QCD calculation). Until recently, s^ has been largely 
unknown and usually assumed to vanish. However, the 
recently published CCFR-NuTeV data on dimuon cross 
sections in i^N and PN scatterin g yi eld a direct handle on 
s{x) and s{x), and hence on s~ because the dimuon 
data reflect semileptonic decays of the charm quark in 



c and W s ^ c events. 



An asymmetric strange sea in the nucleon (s~ 0) 
contributes to a correction term to R~ at LO 01 • If 
the scale dependence of the parton distributions is ne- 
glected, i.e. f{x,Q) ~ f{x), and in the approximation 
of overlooking experimental cuts, the total cross section 
in Eq. (PJ is sensitive to the second Mellin moment inte- 
grals Jdxxfix) of the PDFs HQ- Making the further 
approximation of an isoscalar target, and in the limit of 
a negligible charm quark mass, a strange sea asymmetry 
contributes at LO as 
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where the strangeness asymmetry is quantified by 



[S-] ^ J x[s{x) 



5(x)] dx 
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and [Q^] — J x[q{x)^q{x)]dx withq(2;) — {u{x)+d{x))/2 
represents the isoscalar up and down quark combination. 

By including the dimuon data, and by exploring the 
full allowed parameter space in a global QCD analysis, 
Ref. 13] presents a general picture of the strangeness sec- 
tor of nucleon structure. The strong interplay between 
the existing experimental constraints and the global theo- 
retical constraints, especially the sum rule (|3J), places use- 
ful limits on acceptable values of the strangeness asym- 
metry momentum integral [S~]. The limit quoted in [l^ 
is -0.001 < [S-] < +0.004. A large negative [S'] is 
strongly disfavored by both dimuon and other inclusive 
data. The strict sum rule (jSJ implies that a non-zero 
s~( x) function must change sign at least once. Studies in 
[13| demonstrate that the exact value of [5"] is a volatile 
quantity. The best fit "B" is a solution where negative 
s~ (x) at low X is compensated by positive s~{x) at large 
x; this leads to positivity of the second moment integral 
in Eq. Q . The same trend had previously been observed 
in a fit to inclusive neutrino scattering 4] . Also, this be- 
havior was anticipated by a dynamical model [19| based 
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-0.177 


1.26 
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0.0023 



TABLE I: Shifts in R~ , calculated with PDF sets of 
Ref. [3 (with non-zero [S~]) compared to the value with the 
CTEQ6M set ([5'"] =0), are given in the last column. The 
quality of these new fits is gauged by the relative x'^ values 
(normalized to that of the reference set "B" ) for the dimuon 
data set "2^ and for the subset of global data set which have 
some sensitivity to s~ (x) (labeled "inclusive I"). See )13j for 
details. 

on baryon-meson fluctuations of the nucleon light-cone 
wave function. [27] 

We quantify the impact of the PDFs of Ref. (iSj on the 
Paschos-Wolfenstein relation in Eq. ^ by employing the 
NLO neutrino cross section calculations of Ref. ■ The 
PDF sets A,B,C of Ref. [l^ represent good fits within the 
allowed parton parameter space. They all have s{x) ^ 
s(x), and [5*"] > 0. In our calculations, we employ these 
PDFs consistently; i.e., we use the full set of PDFs, not 
just their strange quark distributions. 

The shift in R~ due to strangeness asymmetry, 5Rj , 
is obtained as the difference: 

SR^ = ^{A3,C,B+,B-} ~ ^CTEQ6 • (^) 

These are given in the last column of Tabled along with a 
summary of the underlying PDFs. We show not only the 
preferred fit values for the sets A, B, C but also results for 
fits that were obtained by using the Lagrange multi- 
plier method to push the limits of the allowed [S~] value 
in both directions somewhat beyond the preferred range 
as described in The quality of the fits is indicated 
by the relative values, which are normalized to the 
reference solution "B". Thus, the values in row B are 1.0 
(italized) by definition. The three preferred sets A,B,C 
are comparable in quality; the extreme sets B"*" and B^ 
are clearly disfavored. 

For a given value of "measured" a shift of the 

theoretical prediction, such as 5Rj , leads to a shift in 
the extracted sin^ 9^ value according to (cf. Eq.QJ): 

5(sin2 Ow) = . (9) 

The results of our calculation (Table 1), along with the 
range -0.001 < [5*-] < 0.004 of Ref. [Taj, which is based 
on more extensive studies than just the fits shown in 
Table 1, lead us to estimate the range of 5R~ , hence 
5{siT? Ow), to be -0.005 < 5{siT? Ow) < -^0.001. 

We find that the shift in R^ , calculated as an aver- 
age over 1/ and V energies according to their fiux spectra, 
is relatively insensitive to the incident neutrino energy. 



The values of 5R- in Table |T| are also approximately 
unchanged when the cut on yEi, is eliminated. These 
findings suggest that the incorporation of other detector 
effects 0, lij , which make the analysis in Ref. [3| some- 
what more involved than a direct measurement of 
will not significantly inroact the importance of the [S~] 
contribution to sin^ 0w |28|. 

The shift in sin^ 9w corresponding to the central fit B 
bridges a substantial part of the original 3 a discrepancy 
between the NuTeV result and the world average of other 
measurements of sin^ 9w- For PDF sets with a shift to- 
ward the negative end, such as —0.004, the discrepancy 
is reduced to less than 1 a. On the other hand, for PDF 
sets with a shift toward the positive end, such as -1-0.001, 
the discrepancy remains. 

More input on s~{x) = s{x) — s{x) would, of course, be 
helpful in pinning down the contribution of strangeness 
asymmetry to 6R~ . Measurements of associated produc- 
tion of charmed jets and bosons at the Tevatron, at 
RHIC, or at the future LHC, would increase our knowl- 
edge of s{x) and s{x) (cf. [131 )• It will help that the "va- 
lence" density s~{x) is more easily accessible than the 
predominantly singlet s{x) + s{x), which is concentrated 
at small x; however, the low expected statistics will make 
this measurement extremely challenging. In principle it 
seems also feasible to study s{x) — s{x) on the lattice [23j |. 
Unfortunately, the most relevant moment [S~] does not 
correspond to a local operator and cannot be calculated 
on the lattice. 

Possible Isospin Violation: Isospin symmetry holds 
to a good approximation in low energy hadron spec- 
troscopy and scattering, but it is not an exact symmetry. 
The level of accuracy of the usual assumption of isospin 
symmetry at the parton level, e.g. Up = d„ and dp = w„, 
is largely unknown. Isospin symmetry violation effects at 
the parton level contribute a shift of the P-W ratio R~ 
by 

""U^e'^" ^^j [Q-] 

where N = {p+n)/2 and, as before, [ ] denotes the second 
Mellin moment. 

There have been model studies jlj that indicate 6RJ 
could be large enough to have an effect on the interpre- 
tation of the NuTeV anomaly. However, it is preferable 
to quantify the allowed range of uncertainty of this effect 
directly and by model-independent global analysis of the 
differences. Unfortunately, there are few experimental 
constraints on these small differences. 

Nonetheless, the MRST collaboration have re- 
cently made a first attempt to separate proton and neu- 
tron PDFs where isospin for the valence quarks is broken 
by a function with a single parameter k. Within phys- 
ically reasonable limits, they find the overall of the 
global fit to be rather insensitive to k. By Eq. H10|l . 
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the determination of sin^ 6w via the measurement of R~ 
is thus subject to a non-negUgible uncertainty due to 
isospin violation. 

To make this point more concrete, we have apphed 
the candidate PDFs from 0| to our NLO calculation, 
in the same spirit as the study of strangeness asymmetry 
discussed above. We find that the range of allowed k 
parameter given in —0.7 < k < 0.7, implies 

- 0.007 < (5i?7 < 0.007 , (11) 

and the best fit value of k = —0.2 corresponds to a shift of 
SRj = —0.0022. A one-parameter functional form may 
not be general enough to pin down the true isospin vi- 
olations of the parton structure. Nevertheless, the large 
range of SRJ in Eq. Ullfl indicates that a reasonable the- 
oretical uncertainty due to isospin violation needs to be 
assigned to the determination of sin^ 

Conclusion: The uncertainties in the parton struc- 
ture of the nucleon that relate to i?" will not decrease 
substantially any time soon. The uncertainties in the 
theory that relates R~ to sin^ 9yj are substantial on the 
scale of precision of the high statistics NuTeV data 0. 
Within their bounds, the results of this study suggest 
that the new dimuon data, the Weinberg angle measure- 
ment, and other global data sets used in QCD parton 
structure analysis can all be consistent within the stan- 
dard model of particle physics. 
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